CYCLOIDEA paralogs function partially redundantly to specify dorsal flower development in Mimulus lewisii.

PREMISE: Duplicated genes (paralogs) are abundant in plant genomes, and their retention may influence the function of genetic programs and contribute to evolutionary novelty. How gene duplication affects genetic modules and what forces contribute to paralog retention are outstanding questions. The CYCLOIDEA(CYC)-dependent flower symmetry program is a model for understanding the evolution of gene duplication, providing multiple examples of paralog partitioning and novelty. However, a novel CYC gene lineage duplication event near the origin of higher core Lamiales (HCL) has received little attention. METHODS: To understand the evolutionary fate of duplicated HCL CYC2 genes, we determined the effects on flower symmetry by suppressing MlCYC2A and MlCYC2B expression using RNA interference (RNAi). We determined the phenotypic effects on flower symmetry in single- and double-silenced backgrounds and coupled our functional analyses with expression surveys of MlCYC2A, MlCYC2B, and a putative downstream RADIALIS (MlRAD5) ortholog. RESULTS: MlCYC2A and MlCYC2B jointly contribute to bilateral flower symmetry. MlCYC2B exhibits a clear dorsal flower identity function and may additionally function in carpel development. MlCYC2A functions in establishing dorsal petal shape. Further, our results suggest an MlCYC2A-MlCYC2B regulatory interaction, which may affect pathway homeostasis. CONCLUSIONS: Our results suggest that CYC paralogs specific to higher core Lamiales may be selectively retained for their joint contribution to bilateral flower symmetry, similar to the independently derived CYC paralogs in the Lamiales model for bilateral flower symmetry research, Antirrhinum majus (snapdragon).

Recent parallel speciation in Antirrhinum involved complex haplotypes and multiple adaptive characters.

A role of ecological adaptation in speciation can be obscured by stochastic processes and differences that species accumulate after genetic isolation. One way to identify adaptive characters and their underlying genes is to study cases of speciation involving parallel adaptations. Recently resolved phylogenies reveal that alpine morphology has evolved in parallel in the genus Antirrhinum (snapdragons): first in an early split of an alpine from a lowland lineage and, more recently, from within the lowland lineage to produce closely related sympatric species with contrasting alpine and lowland forms. Here, we find that two of these later diverged sympatric species are differentiated by only around 2% of nuclear loci. Though showing evidence of recent gene flow, the species remain distinct for a suite of morphological characters typical of earlier-diverged alpine or lowland lineages and their morphologies correlate with features of the local landscape, as expected of ecological adaptations. Morphological differences between the two species involve multiple, unlinked genes so that parental character combinations are readily broken up by recombination in hybrids. We detect little evidence for post-pollination barriers to gene flow or recombination, suggesting that genetic isolation related to ecological adaptation is important in maintaining character combinations and might have contributed to parallel speciation. We also find evidence that genes involved in the earlier alpine-lowland split were reused in parallel evolution of alpine species, consistent with introgressive hybridisation, and speculate that many non-ecological barriers to gene flow might have been purged during the process.

Cloning and functional analysis of RADIALIS-like 1 gene from Antirrhinum majus.

MYB is one of the largest transcription factor families in plants. Among them, the R3-MYB transcription factor RADIALIS (RAD) plays a very important role in the flowers development in Antirrhinum majus. In this study, a R3-MYB gene similar to RAD was found by analyzing the genome of A. majus, which was named AmRADIALIS-like 1 (AmRADL1). The gene function was predicted through bioinformatics. The relative expression levels in different tissues and organs of wild-type A. majus were analyzed by qRT-PCR. AmRADL1 was overexpressed in A. majus, and the transgenic plants were analyzed by morphological observation and histological staining. The results showed that the open reading frame (ORF) of AmRADL1 gene was 306 bp in length, encoding 101 amino acids. It has typical SANT domain, and the C-terminal contains a CREB motif, which was highly homologous to tomato SlFSM1. The results of qRT-PCR showed that AmRADL1 was expressed in roots, stems, leaves and flowers, and the expression level was higher in flowers. Further analysis of its expression in different floral organs showed that AmRADL1 had the highest expression in carpel. The results of histological staining analysis of the transgenic plants showed that compared with the wild type, although the size of the carpel cells of the transgenic plants did not change significantly, the placenta area in the carpel became smaller and the number of cell decreased. In summary, AmRADL1 may be involved in the regulation of carpel development, but the specific mechanism of action in carpel remains to be further studied.

The Snapdragon Genomes Reveal the Evolutionary Dynamics of the S-Locus Supergene.

The genus Antirrhinum has been used as a model to study self-incompatibility extensively. The multi-allelic S-locus, carrying a pistil S-RNase and dozens of S-locus F-box (SLF) genes, underlies the genetic control of self-incompatibility (SI) in Antirrhinum hispanicum. However, there have been limited studies on the genomic organization of the S-locus supergene due to a lack of high-quality genomic data. Here, we present the chromosome-level reference and haplotype-resolved genome assemblies of a self-incompatible A. hispanicum line, AhS7S8. For the first time, 2 complete A. hispanicum S-haplotypes spanning ∼1.2 Mb and containing a total of 32 SLFs were reconstructed, whereas most of the SLFs derived from retroelement-mediated proximal or tandem duplication ∼122 Mya. Back then, the S-RNase gene and incipient SLFs came into linkage to form the pro-type of type-1 S-locus in the common ancestor of eudicots. Furthermore, we detected a pleiotropic cis-transcription factor (TF) associated with regulating the expression of SLFs, and two miRNAs may control the expression of this TF. Interspecific S-locus and intraspecific S-haplotype comparisons revealed the dynamic nature and polymorphism of the S-locus supergene mediated by continuous gene duplication, segmental translocation or loss, and TE-mediated transposition events. Our data provide an excellent resource for future research on the evolutionary studies of the S-RNase-based self-incompatibility system.

Investigation of the cell structure and organelles during autolytic PCD of Antirrhinum majus "Legend White" petals.

One of the classes of the plant developmental programmed cell death (PCD) is vacuolar cell death or autolysis. The results of the transmission electron microscope (TEM) studies indicated that this type of PCD occurs during the petal senescence of Antirrhinum majus "Legend White" flowers. The major hallmarks of the process related to the ultrastructure of the cells involved chloroplast degradation, vacuolation, chromatin condensation, cell wall swelling, degradation of Golgi apparatus, protoplasmic shrinkage, degradation of the endoplasmic reticulum, nuclear fragmentation, rupture of tonoplast, and plasma membrane. Macroautophagy and microautophagy processes were also clearly observed during vacuole formation. As in yeasts, in the present study, Golgi apparatus became autophagosome-like structures during degradation that had autophagy activity and then disappeared. Our results revealed a type of selective microautophagy, piecemeal microautophagy of the nucleus (PMN), in nuclear degradation during PCD of petals that has not previously been reported in plants. Moreover, vesicular structures, such as paramural and multilamellar bodies, were observed in some stages.

How to establish a mutually beneficial relationship between a transposon and its host: lessons from Tam3 in Antirrhinum.

The transposon Tam3 of Antirrhinum (snapdragon) has acquired properties that distinguish it from other transposons. Mobile DNA, commonly referred to as a transposable element or transposon, is considered to be synonymous with a selfish factor. That is, a transposable element increases in copy number and moves copies of itself independently of the survival of the host organism. Therefore, the host collectively regulates the transposition activities of most transposable elements in its genome by epigenetic means. However, our analyses of the structure and behavior of Tam3, as shown by the following five results, provide evidence that it does not behave in a selfish manner in relation to the host. 1) Active transposable elements normally increase the abundance of their non-autonomous elements, whereas Tam3 is known to have no non-autonomous elements, and a limited number of around 10 copies of autonomous elements present in the genome have been isolated as active copies. 2) Tam3 does not transpose at 25 ℃, which is the optimal growth temperature for Antirrhinum. Transposition of Tam3 occurs only at low temperatures of about 15 ℃, which is stressful for Antirrhinum. 3) Few strains of Antirrhinum have been found to contain genes that specifically suppress Tam3 transposition. 4) Most of the Tam3 insertions found in Antirrhinum genes do not affect the host genome, and the expression of these host genes is not completely suppressed. 5) Transcription and translation of the Tam3 transposase gene are not epigenetically regulated by the host. These five experimental results constitute evidence that Tam3 retains features that are dissimilar to those of many other transposons and that it does not behave in a selfish manner that is detrimental to the survival of the host. In this review, we consider what kinds of behavior are required if transposons are to establish a mutually beneficial relationship with their hosts, with reference to Tam3.

Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus.

Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g2. We calculated g2 together with FST and pairwise relatedness (Fij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise Fij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g2. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g2 is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g2, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity.

Digest: Natural selection fluctuates at an extremely fine spatial scale inside a wild population of snapdragon plants.

Marrot et al. used snapdragon plants on a small island to experimentally investigate how spatial structure influences the evolution of biological communities. Using a spline-based fitness function, they studied the varying relationships between traits under selection and driving environmental factors in snapdragons. The authors found that environmental heterogeneity, even on a small spatial scale, may provide several fitness optima on the fitness landscape, paving the way for coexistence of diverse phenotypes. In the absence of sufficient gene flow, this could also lead to microgeographic adaptations.

A phylogeny of Antirrhinum reveals parallel evolution of alpine morphology.

Parallel evolution of similar morphologies in closely related lineages provides insight into the repeatability and predictability of evolution. In the genus Antirrhinum (snapdragons), as in other plants, a suite of morphological characters are associated with adaptation to alpine environments. We tested for parallel trait evolution in Antirrhinum by investigating phylogenetic relationships using restriction-site associated DNA (RAD) sequencing. We then associated phenotypic information to our phylogeny to reconstruct the patterns of morphological evolution and related this to evidence for hybridisation between emergent lineages. Phylogenetic analyses showed that the alpine character syndrome is present in multiple groups, suggesting that Antirrhinum has repeatedly colonised alpine habitats. Dispersal to novel environments happened in the presence of intraspecific and interspecific gene flow. We found support for a model of parallel evolution in Antirrhinum. Hybridisation in natural populations, and a complex genetic architecture underlying the alpine morphology syndrome, support an important role of natural selection in maintaining species divergence in the face of gene flow.

Non-reproducible signals of adaptation to elevation between open and understorey microhabitats in snapdragon plants.

Experimental studies on local adaptation rarely investigate how different environmental variables might modify signals of adaptation or maladaptation. In plant common garden experiments, signals of adaptation or maladaptation to elevation are usually investigated in open habitats under full light. However, most plants inhabit heterogeneous habitats where environmental conditions differ. Understorey microhabitats are common and differ in terms of tree shade, temperature, water availability, microbiota, allelochemicals etc. Germination is a fitness-related trait of major importance for the adaptation of plants to contrasted climate conditions. It is affected by shade in snapdragon plants (Antirrhinum majus) and many other plant species. Here, we tested for the reproducibility of signals extrapolated from germination results between open and understorey microhabitats in two parapatric snapdragon plant subspecies (A. m. striatum and A. m. pseudomajus) characterized by a similar elevation range by using common garden experiments at different elevations. Signals observed under one microhabitat systematically differed in the other. Most scenarios could be inferred, with signals either shifting, appearing or disappearing between different environments. Our findings imply that caution should be taken when extrapolating the evolutionary significance of these types of experimental signals because they are not stable from one local environmental condition to the next. Forecasting the ability of plants to adapt to environmental changes based on common garden and reciprocal transplant experiments must account for the multivariate nature of the environment.

Natural selection fluctuates at an extremely fine spatial scale inside a wild population of snapdragon plants.

Spatial variation in natural selection is expected to shape phenotypic variation of wild populations and drive their evolution. Although evidence of phenotypic divergence across populations experiencing different selection regimes is abundant, investigations of intrapopulation variation in selection pressures remain rare. Fine-grained spatial environmental heterogeneity can be expected to influence selective forces within a wild population and thereby alter its fitness function by producing multiple fitness optima at a fine spatial scale. Here, we tested this hypothesis in a wild population of snapdragon plants living on an extremely small island in southern France (about 7500 m2 ). We estimated the spline-based fitness function linking individuals' fitness and five morphological traits in interaction with three spatially variable ecological drivers. We found that selection acting on several traits varied both in magnitude and direction in response to environmental variables at the scale of a meter. Our findings illustrate how different phenotypes can be selected at different locations within a population in response to environmental variation. Investigating spatial variation in selection within a population, in association with ecological conditions, represents an opportunity to identify putative ecological drivers of selection in the wild.

Bumblebee electric charge stimulates floral volatile emissions in Petunia integrifolia but not in Antirrhinum majus.

The timing of volatile organic compound (VOC) emission by flowering plants often coincides with pollinator foraging activity. Volatile emission is often considered to be paced by environmental variables, such as light intensity, and/or by circadian rhythmicity. The question arises as to what extent pollinators themselves provide information about their presence, in keeping with their long co-evolution with flowering plants. Bumblebees are electrically charged and provide electrical stimulation when visiting plants, as measured via the depolarisation of electric potential in the stem of flowers. Here we test the hypothesis that the electric charge of foraging bumblebees increases the floral volatile emissions of bee pollinated plants. We investigate the change in VOC emissions of two bee-pollinated plants (Petunia integrifolia and Antirrhinum majus) exposed to the electric charge typical of foraging bumblebees. P. integrifolia slightly increases its emissions of a behaviorally and physiologically active compound in response to visits by foraging bumblebees, presenting on average 121 pC of electric charge. We show that for P. integrifolia, strong electrical stimulation (600-700 pC) promotes increased volatile emissions, but this is not found when using weaker electrical charges more representative of flying pollinators (100 pC). Floral volatile emissions of A. majus were not affected by either strong (600-700 pC) or weak electric charges (100 pC). This study opens a new area of research whereby the electrical charge of flying insects may provide information to plants on the presence and phenology of their pollinators. As a form of electroreception, this sensory process would bear adaptive value, enabling plants to better ensure that their attractive chemical messages are released when a potential recipient is present.

Dissipation and transformation of the diamide insecticide cyantraniliprole in ornamental snapdragon (Antirrhinum majus).

Dissipation and transformation of cyantraniliprole, a new diamide class of insecticides, were investigated under greenhouse conditions, using snapdragon (Antirrhinum majus) as the model plant. Dissipation of cyantraniliprole in treated leaves was found to be dependent upon application methods (foliar spray versus soil drench) and doses (high versus low dose), with the parent insecticide being the major residue at various sampling points. A high-dose foliar application resulted in pesticide residue of 6.7-23.8 µg/g foliar fresh weight over 8 weeks of treatments, while in soil drench treatment the residue varied from 0.8 to 1.4 µg/g. However, the residue contents were similar between the two application methods at a low application dose. The transformation pathways of cyantraniliprole were primarily intramolecular rearrangements, with IN-J9Z38 being the major metabolite across treatments. Several other metabolites were also identified, some of which were unique to the application methods. Out of total 26 metabolites tentatively identified in this study, 10 metabolites were unique to foliar application, while six metabolites were unique to soil drench. In addition to plant-mediated biotransformation, photodegradation of the parent compound was identified as a potential mechanism in foliar application.

Discrete bHLH transcription factors play functionally overlapping roles in pigmentation patterning in flowers of Antirrhinum majus.

Floral pigmentation patterning is important for pollinator attraction as well as aesthetic appeal. Patterning of anthocyanin accumulation is frequently associated with variation in activity of the Myb, bHLH and WDR transcription factor complex (MBW) that regulates anthocyanin biosynthesis. Investigation of two classic mutants in Antirrhinum majus, mutabilis and incolorata I, showed they affect a gene encoding a bHLH protein belonging to subclade bHLH-2. The previously characterised gene, Delila, which encodes a bHLH-1 protein, has a bicoloured mutant phenotype, with residual lobe-specific pigmentation conferred by Incolorata I. Both Incolorata I and Delila induce expression of the anthocyanin biosynthetic gene DFR. Rosea 1 (Myb) and WDR1 proteins compete for interaction with Delila, but interact positively to promote Incolorata I activity. Delila positively regulates Incolorata I and WDR1 expression. Hierarchical regulation can explain the bicoloured patterning of delila mutants, through effects on both regulatory gene expression and the activity of promoters of biosynthetic genes like DFR that mediate MBW regulation. bHLH-1 and bHLH-2 proteins contribute to establishing patterns of pigment distribution in A. majus flowers in two ways: through functional redundancy in regulating anthocyanin biosynthetic gene expression, and through differences between the proteins in their ability to regulate genes encoding transcription factors.

Phenotypic Response to Light Versus Shade Associated with DNA Methylation Changes in Snapdragon Plants (Antirrhinum majus).

The phenotypic plasticity of plants in response to change in their light environment, and in particularly, to shade is a schoolbook example of ecologically relevant phenotypic plasticity with evolutionary adaptive implications. Epigenetic variation is known to potentially underlie plant phenotypic plasticity. Yet, little is known about its role in ecologically and evolutionary relevant mechanisms shaping the diversity of plant populations in nature. Here we used a reference-free reduced representation bisulfite sequencing method for non-model organisms (epiGBS) to investigate changes in DNA methylation patterns across the genome in snapdragon plants (Antirrhinum majus L.). We exposed plants to sunlight versus artificially induced shade in four highly inbred lines to exclude genetic confounding effects. Our results showed that phenotypic plasticity in response to light versus shade shaped vegetative traits. They also showed that DNA methylation patterns were modified under light versus shade, with a trend towards global effects over the genome but with large effects found on a restricted portion. We also detected the existence of a correlation between phenotypic and epigenetic variation that neither supported nor rejected its potential role in plasticity. While our findings imply epigenetic changes in response to light versus shade environments in snapdragon plants, whether these changes are directly involved in the phenotypic plastic response of plants remains to be investigated. Our approach contributed to this new finding but illustrates the limits in terms of sample size and statistical power of population epigenetic approaches in non-model organisms. Pushing this boundary will be necessary before the relationship between environmentally induced epigenetic changes and phenotypic plasticity is clarified for ecologically relevant mechanisms with evolutionary implications.

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST -FST comparisons.

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (FST ) and quantitative genetic (QST ) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass-related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.

Radial or Bilateral? The Molecular Basis of Floral Symmetry.

In the plant kingdom, the flower is one of the most relevant evolutionary novelties. Floral symmetry has evolved multiple times from the ancestral condition of radial to bilateral symmetry. During evolution, several transcription factors have been recruited by the different developmental pathways in relation to the increase of plant complexity. The MYB proteins are among the most ancient plant transcription factor families and are implicated in different metabolic and developmental processes. In the model plant Antirrhinum majus, three MYB transcription factors (DIVARICATA, DRIF, and RADIALIS) have a pivotal function in the establishment of floral dorsoventral asymmetry. Here, we present an updated report of the role of the DIV, DRIF, and RAD transcription factors in both eudicots and monocots, pointing out their functional changes during plant evolution. In addition, we discuss the molecular models of the establishment of flower symmetry in different flowering plants.

Shared Mutations in a Novel Glutaredoxin Repressor of Multicellular Trichome Fate Underlie Parallel Evolution of Antirrhinum Species.

Most angiosperms produce trichomes-epidermal hairs that have protective or more specialized roles. Trichomes are multicellular in almost all species and, in the majority, secretory. Despite the importance of multicellular trichomes for plant protection and as a source of high-value products, the mechanisms that control their development are only poorly understood. Here, we investigate the control of multicellular trichome patterns using natural variation within the genus Antirrhinum (snapdragons), which has evolved hairy alpine-adapted species or lowland species with a restricted trichome pattern multiple times in parallel. We find that a single gene, Hairy (H), which is needed to repress trichome fate, underlies variation in trichome patterns between all Antirrhinum species except one. We show that H encodes a novel epidermis-specific glutaredoxin and that the pattern of trichome distribution within individuals reflects the location of H expression. Phylogenetic and functional tests suggest that H gained its trichome-repressing role late in the history of eudicots and that the ancestral Antirrhinum had an active H gene and restricted trichome distribution. Loss of H function was involved in an early divergence of alpine and lowland Antirrhinum lineages, and the alleles underlying this split were later reused in parallel evolution of alpines from lowland ancestors, and vice versa. We also find evidence for an evolutionary reversal from a widespread to restricted trichome distribution involving a suppressor mutation and for a pleiotropic effect of H on plant growth that might constrain the evolution of trichome pattern.

A design principle for floral organ number and arrangement in flowers with bilateral symmetry.

The bilateral symmetry of flowers is a striking morphological achievement during floral evolution, providing high adaptation potential for pollinators. The symmetry can appear when floral organ primordia developmentally initiate. Primordia initiation at the ventral and dorsal sides of the floral bud is differentially regulated by several factors, including external organs of the flower and CYCLOIDEA (CYC) gene homologues, which are expressed asymmetrically on the dorso-ventral axis. It remains unclear how these factors control the diversity in the number and bilateral arrangement of floral organs. Here, we propose a mathematical model demonstrating that the relative strength of the dorsal-to-ventral inhibitions and the size of the floral stem cell region (meristem) determines the number and positions of the sepal and petal primordia. The simulations reproduced the diversity of monocots and eudicots, including snapdragon Antirrhinum majus and its cyc mutant, with respect to organ number, arrangement and initiation patterns, which were dependent on the inhibition strength. These theoretical results suggest that diversity in floral symmetry is primarily regulated by the dorso-ventral inhibitory field and meristem size during developmental evolution.